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Spin pumping enables the vertical transport of pure spin current through Al in a Pd/Al/Ni 80 Fe 20 (Py) trilayer 
film, in which the Py acts as a spin battery. The spin current injected into the Al flows through the Al to reach 
the Pd, resulting in the generation of electromotive forces due to the inverse spin Hall effect in the Pd. The 
electromotive forces decreased with increasing thickness of the Al layer. A simple model based on the theory 
by Tserkovnyak et al., allows an estimation of the spin coherence of the vertical spin transport in the Al of 
61 nm. This comparatively short coherence is attributed to a reduction in spin pumping efficiency because 
of the roughness of the Al/Py interface. 

Since spin-polarized electrons flowing in solids show intriguing characteristics in association with spin 
transfer, spin transport, and spin accumulation 1 , the physical properties of flowing spins are actively studied 
as 'spintronics'. In the field of spintronics, spin pumping 2 " 4 has attracted considerable attention as a method 
of dynamical spin injection, which induces a pure spin current in a nonmagnetic metal 5 . Concurrently, the 
inverse-spin Hall effect (ISHE) in paramagnets, which enables the conversion from a pure spin current to a 
charge current, is frequently used to detect a pure spin current 5,6 . Although the number of studies of the transport 
of a pure spin current (not a spin-wave spin current) in a spacer between a ferromagnet and a paramagnet using 
the dynamical method remains limited 7,8 , previous studies have shown that spin pumping has the potential to 
generate a normal pure spin current in a wide range of materials. The previous demonstrations showed spin 
transport in lateral structures by dynamical spin injection 7,8 . However, it is significant that dynamical spin 
transport is studied in a device with various structures that enables spin transport in various materials, such as 
molecules. 

In this work, we demonstrate successful perpendicular spin transport in the Al layer of a Pd/Al/Ni 80 Fe 2 o(Py) 
trilayer film by the spin pumping at room temperature (RT). The spin transport characteristics are investigated 
using a simple model based on the theory by Tserkovnyak et al 4 . The estimated spin coherence length in the Al was 
comparatively shorter (61 nm) than that observed using a lateral spin valve structure (ca. 300 nm at RT 9,10 ). In 
this paper, we clarify the effect of inhomogeneity of the Al/Py interface on the spin injection efficiency, which 
results in a shorter spin coherence length. 

Result 

Sample structure and inhomogeneity of the Al/Py interface. A schematic structure of the sample used in this 
study is shown in Fig. 1(a). The sample consisted of a Pd/ Al/Py trilayer film with lateral dimensions of 1 mm X 
1.9 mm (/X w). The thicknesses of the Pd (d N2 ) an d Py (d FM ) layers were 15 nm, and that of the Al (d Ni ) was 
varied from 5 to 40 nm. The electrical conductivity of the Py, Al, and Pd were cr FM = 2.05 X 10 6 (Qm)~\ cr N1 = 
9.42 X 10 6 (Qm)~\ and er N2 = 4.08 X 10 6 (Qm)~\ respectively, according to four-terminal measurements. Note 
that we fabricated samples by two different methods: (A) The Pd, Al, and Py layers were formed with a substrate 
cooling system and a deposition rate of 0.25 nm/sec, and (B) the Pd, Al, and Py layers were formed without 
substrate cooling at a deposition rate of 0.04 nm/sec. In the following discussion, these two sample types are 
referred to as sample (A) and sample (B). By using the substrate cooling system, the formation of clusters could be 
suppressed. Additionally, in the case of the vacuum thermal evaporation, the faster the deposition rate, the more 
effectively the formation of the clusters is suppressed, as we demonstrate in the following section. Therefore, by 
comparing the results for samples (A) and (B) we can discuss the effect of the Al/Py interfacial roughness, which 
can have a non-negligible effect on the spin injection efficiency. 
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Figure 1 | Sample structure and equivalent circuit model, (a) A schematic illustration of the Pd(Nl)/Al(N2)/Py(FM) trilayer film. M and H denote the 
precessing magnetization and the external magnetic field, respectively, and the magnetization dynamics of the Py induces spin pumping into the Al. The 
generated pure spin current flows in the Al, reaches the Pd, and is converted to a charge current in the Pd by the ISHE. The dashed arrows shown in the Pd 
layer describe spin scattering by the ISHE. (b) An equivalent circuit model for the Pd/Al/Py trilayer film. 



Dynamical spin injection and inverse spin Hall effect. Whereas 
a pure spin current can influence the magnetization of a 
ferromagnet 11 ' 12 , precessing magnetization loses torque by emitting 
a spin current 2 . In a ferromagnetic resonance (FMR) state, the 
magnetization precession around a static external magnetic field 
enables spin pumping, transferring spin angular momentum from 
the ferromagnet to an adjacent nonmagnetic layer. Therefore, spin 
pumping is a reciprocal process of spin-current-induced magneti- 
zation switching 1112 . The dynamics of the moving magnetization 
M(f) = M + m(f) in the ferromagnet layer are described by the 
Landau-Lifshitz-Gilbert equation as 13 , 



dM(t) 
dt 



a dM(t) 

-yM(f)xHeff+— M(t)x— ^ 
Ms at 



(1) 



where M(f)> M, m(£)> Y> H eff , M s , and a are the magnetization, the 
static component of the magnetization, the dynamical component of 
the magnetization, the gyromagnetic ratio, the effective magnetic 
field, the saturation magnetization, and the dimensionless Gilbert 
damping constant, respectively. Because of the second term on the 
right side of equation (1), the precession spirals down to a static 
magnetization axis by the damping effect. In addition, the damp- 
ing constant becomes large when the emission of spin angular 
momentum from the ferromagnet to a connected nonmagnet 
occurs, i.e., after dynamical spin injection into a nonmagnet 4 . The 
enhancement of the damping effect has been theoretically 4 and 
experimentally 14 studied, and serves as a good index for the 
estimation of the spin current density injected into a nonmagnet 
by spin pumping. In the case of Py, Py has a very small magneto - 
crystalline anisotropy, so the magnetization in the Py was aligned 
along the film plane when the external magnetic field H was applied 
in the film plane. When the H and the magnetization precession 
frequency satisfy the FMR condition 15 , the spins parallel to the 
magnetization precession axis are pumped into the Al layer 2 . 

The ISHE was used to detect the pure spin current 5,6 . The ISHE is a 
process that converts a pure spin current into a charge current, 
resulting in an electric voltage (a schematic image is shown as dashed 
arrows in Fig. 1(a)), and the spin current density is estimated from 
the amplitude of the ISHE voltage when the spin Hall angle is known. 
The charge current J c generated by the ISHE is expressed as 16 , 



Jc 



"ISHE I ^ )Js X r7 « 



(2) 



where J s , $ishe> efi> and a are the spin current density, the spin Hall 
angle, the elementary charge, the Dirac constant, and the spin polar- 
ization vector of the spin current, respectively. By measuring the 
electromotive forces in the Pd under the FMR condition of the Py 
film, we can estimate the spin current density in the Pd layer. 



Before starting the experiments, we measured the spin Hall angle 
of the Pd by the spin pumping method with a Pd/Py bilayer film to 
check the quality of the Pd electrode. The amplitude of the ISHE 
voltage at the Pd and the enhancement of the Gilbert damping con- 
stant were estimated to be 6.0 uV, and 7.3 X 10 7 s" 1 , respectively. The 
spin Hall angle of the Pd was estimated to be 0.017 in our study, 
which agrees with the previously reported value 6 and indicates an 
acceptable sample quality. 

Observation of inverse spin Hall effect. Figures 2(a) and (b) show 
FMR spectra dI(H)/dH and the measured voltages, respectively, for 
the Pd/Al(20 nm)/Py film of sample (A) when an external magnetic 
field H was applied at 9 = 0°, 90°, and 180°, as shown in the right side 
of the Fig. 2(a). A clear electromotive force signal from the Pd layer at 
0 = 0° and 180° was observed around the FMR field, H FMR 
(111 mT). The amplitude of the ISHE voltage Vishe can be 
extracted using the following equation 5 , 

^(H) = Vishe— ~^— 2 + Vahe "^"t^l (3) 



(H-H F 



(h —h FM r) 2 +r 2 



where V(H), T, V ISH e> and Vahe correspond to the output voltage, 
the damping constant, the amplitude of the ISHE, and the amplitude 
of the anomalous Hall effect (AHE), respectively. V ISH e(^— 0°) and 
Vishe(0= 180°) were estimated to be 604 nV and -547 nV, 
respectively, using equation (3). Here, the electromotive force due 
to the ISHE changed sign when the external magnetic field was 
reversed. Therefore, we defined V IS he = (Vi S he(^) — Vi S he($+ 
180°))/2 as an effective voltage for eliminating the heating effect, 
for which the polarity of the voltage does not change the sign. It is 
notable that the amplitude of the ISHE voltage disappeared when the 
external magnetic field was applied perpendicular to the film plane, 
although FMR of the Py was successful (H FMR (6= 90°) was 
measured to be 1210 mT). This confirms the successful detection 
of the ISHE voltage in the Pd film. 

Additional evidence for successful spin injection in the Al by the 
dynamical method was obtained. Figure 3(a) shows the effective 
output voltage V ISH e from sample (A) with 20 nm of Al as a function 
of microwave power. The amplitude of the ISHE voltage is propor- 
tional to the microwave power (see Fig. 3(b)), which is reasonably 
explained by a direct- current- spin-pumping model and also indi- 
cates that the applied microwave power (200 mW) was lower than 
the saturation of the FMR absorption for this system 17 . The linearity 
of the amplitude of the ISHE voltage as a function of the microwave 
power and the inversion of the ISHE voltage polarity with magnet- 
ization reversal indicate successful dynamical spin injection and spin 
transport in the Al. 
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Al-thickness dependence of physical parameters. We quantita- 
tively analyze the transport characteristics of the pure spin current 
in the Al. Figure 4(a) shows the Al-thickness dependence of the 
Vishe value estimated for the Pd/Al/Py films at a microwave 
excitation power of 200 mW. The solid red circles and blue crosses 
show V ISH e of samples (A) and (B), respectively. For sample (A), 
^ishe decreases with increasing d N1 , because spin relaxation in the 
Al decreases the number of spins reaching the Pd film. Moreover, the 
decrease in the total resistance of the sample and the backflow of the 
spin current into the Py with increasing Al film thickness further 
decreased V ISH e- For sample (B), the amplitude of the ISHE voltages 
was small compared with those of sample (A) when the thickness of 
the Al was more than 20 nm. 

The spin current injected into the Pd layer through the Al layer, 
7° S2 , decays along they direction (y = 0 is defined as a point on the Al/ 
Pd interface. See Fig. 1(a).) due to spin relaxation, which we assume 
to be described by 18 , 

sinh[(d N2 -y)/l N2 } t0 , . 

smh{d m /X N2 ) 

where 2 N2 is the spin diffusion length of the Pd. Since the charge 
current density J c (y) can be expressed as Jciy) = Oishe (2e/h) Jsi(y\ the 
average charge current density can be described as 6 , 



</c} = 



(2e\ 




'd m \ 






\2l N2 J 



J° 



S2- 



The total resistance of the trilayer can thus be obtained as, 



V' 



Gpdp + (JNl^Nl + °"N2^A 



-de 



(6) 




o EXP.(0°) 
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•- - FIT. 
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Figure 2 | Observation of the ISHE in sample (A), (a) independence of 
the FMR spectra, d/(H)/dH, measured for the Pd/Al(20 nm)/Py film 
under the microwave excitation of 200 mW. The angle of the external static 
magnetic field was set to be 0, 90, or 180 degrees, as shown at the right side 
of the graph, (b) H dependence of the output voltages for the 
Pd/Al(20 nm)/Py film under the microwave excitation of 200 mW. The 
dashed lines show the theoretical fit obtained using equation (3). 



from an equivalent circuit of the Pd/Al/Py system, as shown in 
Fig. 1(b). Using equations (5) and (6), the spin current density 
injected into the Pd from the Al, /° S2 , is given by, 

a F d F + 0x^x1 + a m d m fh 



T° - 

hi ~ 



V' 



(7) 



The spin current density /° S2 was estimated from V ISH e by setting the 
parameter ^ N2 to 9 nm 19 from equation (7). Figure 4(b) shows the 
Al-thickness dependence of the spin current density, J° S2 . 

Meanwhile, the spin pumping generates a spin current from the 
precessing magnetization M(t), so the direct current component of a 
spin current in the Pd film can be expressed as 6 , 



S2 - 



2tt Jo 



47T 



1 



M(f)x 



dM(t) 



dt 



(8) 



where oo and are the angular frequency of the microwaves (co= 
2nf) and the real part of the effective mixing conductance in the Al 
layer in series with the two interfaces (one with Py and the other with 
Pd). The effective magnetic field H eff can be described as H eff (f) = H 
+ H m (f) + h(f), where H = (0,0,H) is the external magnetic field, 
H m (f) = [0,— 4nm Y (t),0] is the dynamic demagnetization, and h(f) = 
(he i(a \0,0) is the external ac field. m y (t) and h are the y component of 
m(t) and the amplitude of the external ac magnetic field, respectively. 
Under these conditions, the direct-current component of the spin 
current, 7° S2 , at the Pd/Al interface can be rewritten from equations 
(1) and (8) as, 



° 200 mW 
126 mW 
< 50 mW 
o 0.2 mW 




-60-30 0 30 60 
H-H t „ (mT) 

FMR v 7 




A;r 0 100 201 
Microwave power (mW) 

Figure 3 | Microwave power dependence of the electromotive force in the 

Pd. (a) independence of the effective output voltage V ISHE for 
Pd/Al(20 nm)/Py film (sample (A)) under various microwave powers, 
(b) Microwave power dependence of V ISHE . The solid circles are the 
amplitude of the ISHE voltages estimated using equation (3). The solid line 
shows a linear fit to the data. 
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Figure 4 | Al-thickness dependence of physical parameters observed in this study. The solid red circles and the blue crosses shown in (a)-(d) are the 
experimentally obtained values for samples (A) and (B), respectively. The error bars represent the standard error (three samples), (a) Al-thickness 
dependence of the amplitude of the ISHE voltages, (b) Al-thickness dependence of the spin current densities (7° S 2)> which reached the Pd through the Al. 
The dashed line shows the result of the theoretical fitting described in detail in the main text, (c) Al-thickness dependence of the Rrms values of the sample 
surface roughness, (d) Al-thickness dependence of the effective magnetization. The dotted line is a guideline showing M s = 0.065 T. (e) A schematic 
illustration of a cross sectional view of a Py layer deposited on a flat Al layer, (f) A schematic illustration of a cross sectional view of a Py layer deposited on a 
rough Al layer. 



4= 



4-KM s y + J (4nM s f)> 2 + 4ct> 2 



87ia 2 [(47tM s ) 2 y 2 +4co 2 ] 



(9) 



In this study, some of the pumped spins are scattered in the Al layer, 
some reach the Pd layer and are converted to a charge current, and 
the remainder diffuses back into the Py layer, which is different to the 
simple bilayer case. Here, the spins diffusing back into the Py layer 
return the spin angular momentum to the magnetization of the Py 
layer. Therefore, because of the backflow, the enhancement of the 
dimensionless Gilbert damping constant is reduced relative to the 
case of the Pd/Py bilayer 4 , which should be taken into account in our 
analysis. The enhancement of the dimensionless Gilbert damping 
constant due to spin relaxation in the trilayer system is given in 
MKS units as 4 , 



<dm) = Ws\ G0 + 



1 ^&N1/FM\ 



1 + tanh(d N i /k m )g^ 2/m \J 3T SF /T e/ (n/ 2k 2 F ) 



(10) 



tanh(d m /X m )+gl[ /m y / 3Ts F /?ei(n/2k 2 F ) 



4nh d F 



where G 0 ,^ 1/FM ,&^ /N p t S f> ?el> /c f , V, and# L , are the bulk Gilbert 
damping constant of Py (G 0 = yM s a 0 ), the real part of the mixing 
conductance of the Al/Py interface, the real part of the mixing con- 
ductance of the Pd/Al interface, the spin relaxation time, the elastic 
scattering time, the Fermi wavenumber, the spin diffusion length of 



the Al, and the g- factor of electrons, respectively. Setting g^, 



2/Nl 



= 0 



cc(d m ) — cc(d m ) n 



(ii) 



Gilbert damping constant of the Al/Py bilayer. Here, the real part of 
the effective mixing conductance, g^, for the Pd/ Al/Py film can be 
determined from the enhancement of the dimensionless Gilbert 
damping constant due to the missing of the spin angular momentum 
during the spin pumping as 20 , 

II 4nM s d F 

Since the difference between the dimensionless Gilbert damping 
constants of Pd/ Al/Py and Pd/Al depends on the presence of the 
Al layer, the effective mixing conductance also depends on the thick- 
ness of the Al layer. Using the amplitude of the observed ISHE 
voltages and equation (7), we obtain the spin current density injected 
into the Pd from the Al. Then, using equations (9), (10), and (11), the 
spin diffusion length of Al, /L N1 ; the mixing conductances, glj[/ FM and 
&N2/NV anc ^ me Gilbert damping constant of the bulk Py film, G 0 can 
be estimated. 

Using equations (9)-(ll) with M s = 0.065 T, k F = 1.75 X 
10 10 m~\ and h = 6.0 X 10" 5 T and a spin flip probability 8 = T d / 
t sf = 1.5 X 10" 3 21 ' 22 , we found G 0 = 7.8 X 10 7 s-\gl\ /FM = 4.0 X 
10 20 m" 2 ,^ N1 = 4.0 X 10 19 m" 2 , and2 N1 = 61 nm for the Pd/Al/ 
Py film (the fitting line is shown as a dashed line in Fig. 4(b)). Here, 
the sample data (Al thickness = 15 nm) is not included for the 
fitting, since its spin current density was singularly large and the 
fitting calculation was not successfully completed, for reasons cur- 
rently under investigation (we have confirmed that all of the experi- 
mental results are reproducible). The Gilbert damping constant of 
the bulk Py film, G 0 , and the real part of the mixing conductance, 
&vi/FM' are c l° se t0 t ^ ie vames reported in the literature 23,24 , and the 
bulk damping constant of the Py film was estimated to be G 0 = 7 X 
10 7 s" 1 on the basis of FMR measurements, and the real part of the 
mixing conductance of Al/Py was reported to be g^/ PAf = 2.0 X 



decouples the Pd/Al junction, i.e., we can obtain the dimensionless 10 20 m 2 . Although the real part of the mixing conductance of the 
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Pd/Al has not yet been reported, the real part of the mixing conduc- 
tance of Cu/Pt was reported to be g^ 2 / N1 — 3.5 X 10 19 m" 24 , which is 
comparable to the value for our Pd/Al. 

Discussion 

The estimated spin diffusion length of the Al in this study is smaller 
than that estimated using an electrical method 910 ' 21,25 . We assumed 
that this was due to the effect of interfacial roughness of the Al/Py 
interface on the spin injection efficiency, the following experiments 
were implemented. Figure 4(c) shows the Al- thickness dependence 
of the root mean square roughness R RM s values of the surface rough- 
ness of samples (A) and (B). The Rrms values of sample (A) are 
smaller than those of sample (B). This difference is especially obvious 
for thicker Al layers (d N1 > 25 nm), which is ascribable to Al cluster 
formation induced by the sample temperature increase by heat radi- 
ation during the evaporation. It is notable that a dramatic decrease in 
spin current density was observed for samples with the same Al 
thickness (d N1 > 25 nm, see Fig. 4(b)). This finding indicates that 
interfacial roughness between the Py and the Al impedes efficient 
spin injection. The effective magnetization, which was estimated 
from the FMR measurements, decreased with increasing Al thickness 
(see Fig. 4(d)), and this result strongly supports the above discussion 
and conclusions. Figures 4(e) and (f) show schematic illustrations of 
cross sectional views of a Py layers deposited on Al films of different 
roughness. For a sufficiently flat Al layer surface, such as that shown 
in Fig. 4(e), the Py film is uniformly magnetized along the in-plane 
direction. Meanwhile, for a rough Al layer surface, such as that 
shown in Fig. 4(f), the magnetization of the Py film is easily tilted 
from the in-plane direction 26 , which decreases the effective magnet- 
ization. The reduction of the effective magnetization is attributed to 
the roughness, and the rough interface impedes efficient spin injec- 
tion through the interface because the spin current density is reduced 
due to the low effective magnetization. Hence, the spin injection 
efficiency decreases with increasing roughness of the Al surface. 
Our experimental results regarding the effective magnetization and 
roughness are in good agreement with the above interpretation. The 
findings in this study indicate that the control of the interface is of 
great significance in achieving effective dynamical spin injection. We 
expect that the combination of spin pumping and the ISHE can be 
used for spin injection, transport, and detection in various systems 
with vertical structures, such as molecular spin valves. 

In conclusion, the Al thickness dependence of the amplitude of the 
ISHE voltage and the spin current density were studied using Pd/Al/ 
Py trilayer films. The Gilbert damping constant of the bulk Py, the 
real part of the mixing conductance of the Al/Py, the real part of the 
mixing conductance of the Pd/Al, and the spin diffusion length in the 
Al were estimated to be G 0 = 7.8 X 10 7 s~\gl[ /FM = 4.0 X 10 20 m" 2 , 
&V2/N1 = x 10 19 ni" 2 , and/l N1 = 61 nm, respectively, by extend- 
ing the conventional theoretical model. This estimation revealed the 
importance of interfacial control in efficient spin injection. The mea- 
surement and estimation methods provide a valuable technique for 
estimating the spin coherence of materials, where successful and 
reproducible spin injection and transport are awaited. 

Methods 

Experimental procedure. Two Au pads were formed as electrodes at both ends of the 
Pd (see Fig. 1(a)) by electron beam deposition using a shadow mask method. The Pd 
film, Al film, and Py film were formed on a thermally oxidized Si substrate by electron 
beam deposition using a shadow mask method with vacuum in-situ process. For 
sample (A), the temperature was kept at — 2°C during the deposition by a cooling 
medium flowing through the sample holder. Two contacts for taking the inverse spin- 
Hall effect (ISHE) voltage measurements were attached to the Au Pads with Ag paste. 
The electromotive force due to the ISHE is obtained as Vishe = [RfRniRni/ '(RpR-m + 
#ni#n2 + #n2#f)Uc as shown in Fig. 1(b). I c = /d N2 </c> is charge current 
generated by the ISHE. In the FMR measurements, each sample was placed at the 
center of a TE 102 microwave cavity in an electron spin resonance system (BRUKER, 
EMX10/12 EPR spectrometer) with a frequency off = 9.62 GHz at room 
temperature. In this system, the microwave amplitude h was 6.0 X 10" 2 mT for a 
microwave power P M w of 200 mW. An external magnetic field H was applied at 



6 = 0°, 90°, and 180°, as shown in Fig. 1(a). We used the resonance condition, 
{(v/y) 2 = H F mr(H FM r + 47iM s ). 
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